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ABSTRACT 

Context. We study the atmosphere of the carbon-rich Mira RU Vir using the mid-infrared high spatial reso¬ 
lution interferometric observations from VLTI/MIDI. 

Aims. The aim of this work is to analyse the atmosphere of the carbon-rich Mira RU Vir, with hydrostatic 
and dynamic models, in this way deepening the knowledge of the dynamic processes at work in carbon-rich 
Miras. 

Methods. We compare spectro-photometric and interferometric measurements of this carbon-rich Mira 
AGB star, with the predictions of different kinds of modelling approaches (hydrostatic model atmospheres 
plus MOD-More Of Dusty, self-consistent dynamic model atmospheres). A geometric model fitting tool is 
used for a first interpretation of the interferometric data. 

Results. The results show that a joint use of different kind of observations (photometry, spectroscopy, inter¬ 
ferometry) is essential to shed light on the structure of the atmosphere of a carbon-rich Mira. The dynamic 
model atmospheres fit well the ISO spectrum in the wavelength range A = [2.9,25.0] pm. Nevertheless, a 
discrepancy is noticeable both in the SED (visible), and in the interferometric visibilities (shape and level). 
A possible explanation are intra-/inter-cycle variations in the dynamic model atmospheres as well as in the 
observations. The presence of a companion star and/or a disk or a decrease of mass loss within the last few 
hundred years cannot be excluded but these explanations are considered unlikely. 

Key words, instrumentation: high angular resolution - techniques: interferometric - stars: AGB and post- 
AGB - stars: atmospheres - stars: circumstellar matter - stars: fundamental parameters 


1. Introduction 

The Asymptotic Giant Branch (AGB) is the late 
evolutionary stage of low-to-intermediate mass stars 
(typically ss 0.8 to ss SMq). These objects are char¬ 
acterized by a C-O core and He/H burning shells, 
surrounded by a convective envelope, while the at- 


* Based on observations made with ESO telescopes at 
La Silla Paranal Observatory under program IDs: 085.D- 
0756 and 093.D-0708 


mosphere consists of atomic and molecular gas, in¬ 
cluding dust grains. Stellar pulsation can generate 
shock waves that run through the atmosphere. These 
shocks can propagate outward, causing a levitation 
of the outer atmosphere layers, improving the condi¬ 
tions for the formation of dust, which may then lead 
to a dust-driven wind. 

At the early stages of their life, AGB stars have a 
carbon-to-oxygen-ratio below one. The third dredge- 
up can turn the chemistry of these objects from 
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oxygen- into carbon-rich (Iben & Renzini 1983| l. 
Carbon-rich AGB stars are important contributors 
to the enrichment of the interstellar medium, and 
in their spectra there are signs of carbon-bearing 
molecules such as C 2 , C 3 , C 2 H 2 , CN or HCN, while 
the dust is mainly dominated by amorphous carbon 
dust grains and SiC ( Loidl et al.|2001[|Yamamura & 
|de Jong|2000| l. 

By studying the stellar atmospheres of AGB 
stars, we hope to better understand many processes 
that can occur, e.g. the dust formation and mass 
loss via strong stellar winds, the connection be¬ 
tween pulsation and atmospheric structure. In the 
case of C-rich AGB stars with no pronounced pul¬ 
sation, most of the observables derived from hydro¬ 
static models agree fairly well with the observations 
( Aringer et al. 12009] ). However, as the star evolves, 
time-dependent processes become more important, 
and the atmosphere expands, causing a decrease of 
the effective temperature. Hence, model atmospheres 
that take into account those dynamic processes (pul¬ 
sation, dust formation, mass loss) are necessary (e.g. 


1997t IHbfner et al.||2003[ see also the excellent re 

views from Woitke|2003 and |Hofner|2007| l . 

Over the last years there has been an increase 
in the number of works that investigate the atmo¬ 
sphere of AGB stars combining various techniques 
dWittkowski et al.||20m1 |2008l |20TT1 |Neilson &| 
|Lester||2008[ [Marti- Vidal et al.||2011| l. On the other 
hand, to date only few interferometric observations 
of carbon stars have been compared with model at- 
mospheres (jSacuto et al.|2011| Paladini et al.|2011 


Cruzalebes et al.||20131 |Klotz et al.|2013t|van Belle| 


etal.|2013D . 

In this paper we intend to study the atmosphere 
of the carbon-rich Mira RU Vir, by means of pho¬ 
tometry, spectroscopy, interferometry, and model at¬ 
mospheres (jMattsson et al.[ 2010[ Aringer et al. 


20091, using the most recent grid of dynamic atmo¬ 


sphere models and synthetic spectra (Eriksson et al. 

RU Vir is a C-rich AGB Variable Star of Mira 
type. It has a period of 433.2 days in the V band 
(Samus et al.|2009]). Its distance is 910 pc (based on 


the Period-Luminosity relation from Whitelock et al. 


2006 1 , the amplitude of the variability in the V band 


is 5.2 mag and the average value of the magnitude 
in y is < y >= 11.6 mag. RU Vir is surrounded by 
a carbon-rich dusty shell made of amorphous carbon 
(AmC) and silicon carbide (SiC). The latter shows its 


Bowen|19881|Fleischer et al.|1992[|Hofner & Dorfi| 


presence through the strong emission feature around 

~ 11 jum. 

We present MIDI (MID-infrared Interferometric 
instrument) observations (P.I. G. Rau: ID 093.D- 
0708(A)) probing regions of dust formation in the 
circumstellar envelope ( pvfowotny et al. 2011 Danchi 
et al. |1994| l. In a forthcoming paper, we will also 
consider a larger sample of carbon-rich AGB stars. 

In Sect. 1^ we present the photometry, the spec¬ 
troscopy, and the interferometric observations of 
RU Vir and describe the data reduction performed. 
Sectionj^explores the geometry of the environment, 
and the fit of the geometrical models to the interfer¬ 
ometric visibilities to constrain the morphology and 
the brightness distribution of the object. In general 
the visibility V as a function of spatial frequencies 
u and V, as described e.g. in |Winters et al.|p995| l, is 
the two-dimentional Fourier transform of the inten¬ 
sity distribution /(t;, ^), where 77 and ^ are the corre¬ 
sponding coordinates on the celestial sphere. In the 
special case of a circular symmetric intensity distri¬ 
bution, which is produced by spherical symmetric 
models, the visibility is given by the Hankel trans¬ 
form of the intensity profile: 


y{q) 


r<l> 

Jo 


I(<p) ^Jo{27T<pq) dcf) , 


where q = is the spatial frequency, (f) is 

the angular separation from the center of the star 
(0 = rjd) with d the distance of the star, r the radius 
and Jq the zeroth-order Bessel function. Therefore, 
the term interferometric visibilities will be from now 
on referred to as visibilities V, used to denote the nor¬ 
malized quantity V{q)IV{Q). In Sect, ^an d|5] we in¬ 
vestigate the modelling of the atmosphere of RU Vir, 
from the hydrostatic and hydrodynamic point of 
view. In Sect.|^we propose some of the possible sce¬ 
narios that could justify the optical excess (Sect. 6.2 1 


and the discrepancy in the visibilities (Sect. 6.3 1 . 


Moreover, a discussion on the stellar parameters can 
be found in Sect. 6.1 Lastly, we present the conclu¬ 
sions and perspectives for future works, in Sect. [7] 

2. Observationoal data for RU Vir 

2.1. ISO spectrum 

RU Vir was observed with the Short Wavelength 
Spectrometer on board the ISO satellite (SWS, de 
Graauw et al.jl996] l once on 20 July 20, 1996 (Sloan 
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et al.| 20031. This ISO/SWS spectrum covers the 


wavelength range form 2.36 to 45.35 fim, of a spec¬ 
tral resolution of R a; 200. To the purpose of our in¬ 
vestigations, an error of + 10% is assumed for wave- 
lengt < 4.05 /tm, and + 5% towards redden wave¬ 
lengths (Sloan et al. 2003 1 . The ISO spectrum can 
be seen in black in Fig.|^ 


2.2. Photometry 

RU Vir photometry in the filters B, V, R, I, J, H, K, is 


|et al. 

(2011 

1 , we pick the B, V, R, I observations 

from 

Eggen 

(|1975|l and J, H, K observations from 

Whitelock et al.|(|2006i. The photometry is taken at 


the same phase (p as the ISO/SWS spectrum, which 
can be calculated in the following way: 


j,RU Vir 
USO 


(f-Wl 


= 0.65 + 0.1 


( 1 ) 


where t is the time of the observation, in Julian 
date. To is the selected phase-zero point, correspon- 
ing to the nearest visual light maximum of the star. 
The light curves were interpolated to this phase 
( [Nowotny et al.|[T011[ Fig. 8 there). The uncertain¬ 
ties of the different photometric measurements are 
estimated from the uncertainty in the phase and the 
scatter of the phased light curve data. 

The visual light curve of RU Vir from AAVSO 
observations is shown in Fig|T] The time of the ob¬ 
servation of the ISO spectrum and the time interval 
of |Eggen ( 1975| l photometry (that covers more than 
one pulsation cycle) are underlined in vertical green 
and blue lines, respectively. The light curve favor 
the presence of a secondary period on a timescale 
of ~ 30yr ( Percy & Bagby|19^ . 


2.3. MIDI data 

RU Vir was observed in 2010 (ID: 085.D-0756(C)) 
and 2014 (ID: 093.D-0708(A)) with the 1.8m 
Auxiliary Telescopes (ATs) of the Very Large 
Telescope Interferometer MIDI ( Leinert et al.|2003) l. 
MIDI provides wavelength dependent visibilities, 
photometry, and differential phases in the V-band 
(drange = [8, 13] jjm). 

The journal of the observations is available in 
Table A.l of the appendix. The calibrators used 
are listed below the corresponding science, and the 
second-last column identifies if the observations are 


Table 2: Parameters of the calibrator targets. 


Target 

Spectral type 

[Jy] 

Diameter 

[mas] 

HD120323 

M4.5m 

255.4 

13.25 ± 0.060 

HD133216 

M3/M4III 

200.7 

11.154 + 0.046 

HD81797 

K3II-III 

157.6 

9.142 + 0.045 

Notes. (a)IRAS Point 

Source 

Catalog: http: 

//simbad.u- 

strasbg.fr/simbad/ 


(b] http: //WWW. ster. kuleuven. ac 

.be/~tijl/MIDI_ 


calibration 


carried out in SCI-PHOT or HIGH-SENS mode. All 
the observations are in low spectral resolution (R = 
30). The Mv-coverage is shown in Eig. The main 
characteristics of the calibrators are listed in Tablej^ 


The choice of the calibrators follows the selec¬ 
tion criteria described in |Klotz et al.|P012al l. 

Data were reduced with the latest data reduction 
software package MIA+EWS ( Jaffe|20()4) l. The cali¬ 
brated visibilities at different baselines are plotted in 
Eig. The wavelength dependent visibility exhibits 
the typical shape of the visibilities for carbon stars 
with dust shells containing SiC grains (Paladini et 
al. in prep.). There is a drop between 8-9 fim due to 
C 2 H 2 + HCN opacities, and another decline in the 
visibility shape at ~ 11.3 /urn, due to SiC. 


The difference in position angle between the two 
sets of the observations (in 2010 and 2014) does not 
allow us to perform a check for interferometric vari¬ 
ability because they do not probe the same spatial 
frequencies. So far, interferometric variability was 
observed only for one star (V Oph, Ohnaka et al.| 
2007|l; therefore we decided to analyze the data all 


at the same time. When it comes to study the geom¬ 
etry (Sect. ^ we will keep in mind that if these two 
points push for an asymmetric solution, then vari¬ 
ability could be still an explanation. 


The MIDI differential phase is zero, and the 
MIDI spectra are shown in Eig. The shape of the 
MIDI spectra agree within the errors with the ISO 
spectrum. The difference in flux level between the 
two sets of MIDI observations is typical for C-rich 
Miras ( |Hron et al.|1997| ). 
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Table 1: RU Vir photometry in different filters adopted from Nowotny et al. (2011 1 , in units of mag. The 


photometry is interpolated to the phase of the ISO/SWS spectrum (see Sect. 2.21. 


B 

V 

R 

/ 

J 

H 

K 

15.95 ± 0.90 

11.75 + 0.70 

8.50 + 0.50 

6.90 + 0.50 

5.50 + 0.20 

3.75 + 0.15 

2.30 + 0.20 





BVRI JHK 

MIDI 
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Fig. 1: AAVSO light curve of RU Vir in black. The gray shadow shows the time range of the BVRI pho¬ 
tometry observed by |Eggen] ( | 19751 1. The vertical green, blue, and violet lines denote the epoch of the ISO 
spectrum, JHK photometry, MIDI observation respectively. 


3. Geometrical models 

3. 1. Geometry of the environment 

Before comparing the model atmospheres and the 
data, we studied the morphology of the circumstellar 


environment, interpreting the MIDI interferometric 
data with geometric models. 


This has been done by using the geometrical 
model fitting tool GEM-FIND (GEometrical Model 


Eitting for INterferometric Data) of Klotz et al. 


4 






















































Rau et al.: Modelling the atmosphere of the carbon-rich Mira RU Vir 
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Fig. 2: Mv-coverage of the MIDI RU Vir observations 
listed in Table 


A.l 


dispersed in wavelengths. 



Wavelength [/vm] 


Fig. 3: RU Vir calibrated visibilities, at different 
baselines and projected angles. 


|2012b| The program fits geometrical models to 
wavelength dependent visibilities in the V-band, to 
constrain the morphology and brightness distribu¬ 
tion of an object. A detailed description of the fitting 
strategy and of the minimization procedure can 
be found in |Klotz et al.| ( |2012al l. 

Different models with optically thick (models 
with one-component) or thin (models with two- 
components) circumstellar envelopes were tested. In 
Tablej^we list the corresponding reducedvalues, 
and the parameters fitting best. In the case of elliptic 
models the ratio of major and minor axis is shown. 


Fig. 4: Comparison of observational data with mod¬ 
elling results. The black line shows the ISO/SWS 
spectrum of RU Vir, the green line the SED corrected 
from dust emission. The orange and violet shadows 
marks the MIDI spectra, while the gray one denote 
hydrostatic models grid. The best fitting model is 
shown in cyan. 


For the two component models with circular uniform 
disc (UD) -H circular gaussian, the brightness ratio 
is indicated. In the latter case the fit was performed 
with a fixed value of the angular diameter, to avoid to 
have too many free parameters in the fit. We started 
the fit of the composite model imposing the diameter 
to be equal to the one derived from the 0/(V - K) re¬ 
lation of |van Belle et al.|P013| l: 9 - 2.18 mas. Then 
we increased the diameter by steps of 2 mas until a 
f'max =20 mas. 

One component spherically symmetric and ellip¬ 
tical models have large reduced and are not able 
to reproduce the data. A better fit to the data is ob¬ 
tained by a two-components model, including a cir¬ 
cular UD plus a gaussian distribution of the circular 
dusty envelope. This yields a reduced x^i^ 
for a central star simulated by a UD with angular di¬ 
ameter at 11 fim of 18 mas (see Table |^, probably 
because of molecular opacities. 

In Fig. we show the the visibilities vs. wave¬ 
length, while in Fig.j^the fit of the visibility vs. base¬ 
line at different wavelength is shown. 

In general, the geometric models show that, for 
the u-v covered by our RU Vir observations, there is 
no major asymmetry detected. 


5 






















































Rau et al.: Modelling the atmosphere of the carbon-rich Mira RU Vir 

Table 3: GEM-FIND results. For each fit (one or two components), the full width at half maximum and an¬ 
gular diameter at 11 /rm are given. When applicable also the ratio of minor and major axis and the brightness 
ratio are given. The best fitting model is shown in bold face. 


Geometric model 

^red 

FWHMu 

[mas] 

dn 

[mas] 

a/b 

B ratio 

One component 

Circular UD 

40.10 


38.10 ±0.17 



Circular Gauss 

15.42 

28.44 ± 0.24 




Elliptic UD 

12.31 


177.34 ±0.77 

0.2 


Elliptic Gauss 

5.41 


56.62 ± 0.44 

0.3 


Two components 

Circ UD-l-Circ Gauss 

1.59 

37.88 ±0.99 

2.81 “ 


0.12 ±0.01 

Circ UD-l-Circ Gauss 

1.57 

38.02 ± 1.01 

6.00“ 


0.12 ±0.01 

Circ UD-tCirc Gauss 

1.53 

38.18 ± 1.02 

8.00“ 


0.13 ±0.01 

Circ UD-tCirc Gauss 

1.49 

38.38 ± 1.05 

10.00 “ 


0.13 ±0.01 

Circ UD-tCirc Gauss 

1.43 

38.65 ± 1.08 

12.00 “ 


0.15 ±0.01 

Circ UD-l-Circ Gauss 

1.35 

39.01 ± 1.12 

14.00 “ 


0.15 ±0.01 

Circ UD-l-Circ Gauss 

1.27 

39.47 ± 1.18 

16.00 “ 


0.17 ±0.01 

Circ UD-l-Circ Gauss 

1.15 

40.08 ± 1.12 

18.00“ 


0.19 ±0.01 

Circ UD-tCirc Gauss 

1.01 

40.90 ± 1.38 

20.00 “ 


0.21 ±0.01 


Notes, (a): Models with central star diameter fixed at the values indicated above during the fit. 


X = 8.5 [/Urri] 



X = 1 1.5 [^irn] 



model atmospheres and 


described in Aringer et 
of 


the computed spectra are 
. The resolution 
R 


runrjTsR 


Fig. 5: Comparison of observational data with the 
best fitting geometrical model. The calibrated vis¬ 
ibilities are shown versus baseline length, for two 
different wavelengths (left; 8.5 jum; right; 11.5 fim). 
The symbols represent MIDI observations at three 
different baseline configurations. The lines show the 
best-fitting model consisting of a circular UD and a 
circular gaussian. 


4. Hydrostatic model atmospheres + MOD 
models of dusty envelopes 


Following Sacuto et al.| ( |201l] l, we perform the first 
step of our study first by fitting the ISO spectrum 
with hydrostatic model atmospheres (COMARCS 
models, [Aringer et al.|2009 1. 

Each hydrostatic model is identified by the fol¬ 
lowing initial parameters; mass, effective temper¬ 
ature, C/O, log(g), metallicity. Details about these 


each synthetic spectrurrlM is R = 200 to 

match the resolution of ISO, and the spectral range 
[0.4,45.0] ^lm. 

Since RU Vir is a dust-enshrouded carbon-Mira, 
we included dust by using the radiative transfer code 
MOD (More Of Dusty, Groenewegen 2012| ). This 
code is based on DUSTY, that is a publicly available 
1 D dust radiative transfer code ( jlvezic & Elitzurj 
|1995| l, whose mathematical formulation is described 
in detail in jlvezic et al.j ( [T997| l. MOD tries to opti¬ 
mize a set of parameters (e.g. luminosity, dust opti¬ 
cal depth, dust condensation temperature and slope 
of the density law) with some constraints as spectro- 
photometric-data, ID intensity profiles and visibil¬ 
ity curves. A quantitative measure of the quality of 
the fit is obtained through performing a;^^^ minimiza¬ 
tion, where the;^^^ are computed separately for all the 
four types of observations. We note at this point that 
the dust condensation temperature in this context ac¬ 
tually is the dust temperature at the inner bound¬ 
ary of the dust shell. Following previous work (e.g. 
Groenewegen|20 1 2)l we will call it the dust conden¬ 


sation temperature. 

' Available at 
synphot/Cstars/ 


http://stev.oapd.inaf.it/ 
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Following [Sacuto et al.| (2011l, we use a mix¬ 
ture of 90 % AmC and 10 % SiC for the dust 
composition. This ratio gave the best results for the 
ISO/SWS spectrum of the carbon star R Scl. Fits for 
RU Vir with a lower SiC fraction were not as good. 
The optical constants for opacities of silicon carbide 


and amorphous carbon are taken from Pitman et al. 
( [2008| l and [Rouleau & Martin] ( |1991| ) respectively. 
As the standard grain size distribution of DUSTY 
produced much too peaked SiC features, we used the 
Distribution of Hollow Spheres (DHS, jOroenewegen] 


20 12] and references therein) with a mean grain size 


of the dust mixture of 0.2 /rm. This grain size is 
similar to the typical sizes found from models for 
dust driven mass loss ( [Matts son & H6fner||201 ij l. 
However, we note that these models include only 
amorphous carbon grains. On the other hand, SiC 
grains found in presolar meteorites have sizes of the 
order of 1 jum (e.g. Gail & Sedlmayr|2013| l, indicat¬ 
ing that our mean size is consistent with laboratory 
data. The chosen grain size distribution then, gave a 
good fit of the SiC feature. 


4.1. The fitting procedure 


The hydrostatic models alone are not able to repro¬ 
duce the ISO spectrum. This is because the hydro¬ 
static models do not include dust, while the dust con¬ 
tent in the atmosphere of RU Vir is quite prominent, 
as the spectral energy distribution from ISO/SWS 
shows (see Fig. black line). In the following we 
will use the hydrostatic model atmosphere as inner 
radiation source for MOD. The model parameters 
are selected as described below. Considering that 
RU Vir is located in the vicinity of the sun, a first 
selection of the hydrostatic models was made choos¬ 
ing only models with solar metallicity. To the aim of 
being independent of the distance each hydrostatic 
model spectrum was normalized to the ISO flux at 
2.9 pm. We chose this wavelength because here re¬ 
sides the local minimum of the molecular absorption 
( [Aringer et al.|20()9] l. 


We constrained the effective temperature of 
RU Vir by fitting the ISO spectrum with the syn¬ 
thetic spectra, and performing a minimization, in 
the range [2.9,3.6] pm. The molecular features in 
this wavelength range are good temperature indica¬ 
tors for hydrostatic C-stars pprgensen et al. 2000| l 
and provide a first approximation for our case. The 
result obtained in this way is T = 3100 + 100 K. The 


40 


30 

X 

20 


10 

0 

2000 2500 3000 3500 4000 4500 

T,„ [K] 

Fig. 6; Results of the fitting for the effective 
temperature, of the the synthetic spectra based on 
COMARCS models, with the ISO spectrum. The fit 
is done in the range of: [2.9,3.6] ;um. 



C/O is set to C/O- 2.0 because this value is the one 
corresponding to the lowest in the wavelength re¬ 
gion above. Figure [^illustrates the;^^^ vs. in the 
wavelength range chosen as the best for the parame¬ 
ter determination; the best fitting value is located at 
the minimum of the^^. 

We proceed including the dust component with 
MOD a posteriori, in the way described as follows. 
We fit with MOD the photometry (see Sect. 2.2 1 , 
spectroscopy (see Sect. 113 and interferometry (see 
Sect. 2.3 I of the star. The input parameters of MOD 
are effective temperature (determined with the hy¬ 
drostatic model atmosphere fit), distance and the 
dust properties. The optical depth r, the parameter 
that governs the shape of the density law p, and the 
dust condensation temperature T^ can be fitted or set 
to a fixed value. 

To assess how much the spectrum used for the in¬ 
ner radiation source is affected by dust emission, an 
iterative procedure was used. We determined a ’’cor¬ 
rected” ISO spectrum by performing a MOD fit us¬ 
ing the 3100K spectrum (see Fig. [^ derived above. 
We then subtracted the contribution of the dust emis¬ 
sion at each wavelength from the total flux, obtain¬ 
ing in this way a stellar photospheric contribution. 
Repeating the effective temperature determination, 
results in a best model atmosphere with an effec¬ 
tive temperature of 3000 K, i.e. only slightly cooler 
than the first estimate. This model spectrum was then 
used for the remaining calculations. We note that the 
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Table 4; Output values of MOD fitting, for three different cases. The parameters are; luminosity L, optical 
depth T, dust condensation temperature T^, parameter of the slope of the density law p. The last four lines 
show the reduced for the fit of; photometry, spectroscopy, interferometric visibilities, and then the total 
one. 



T; = 1000 K 

T* = 1300 K 

all free 

L[Lo] 

6118.05 ±32.54 

5770 ± 28.02 

6322.66 ±51.07 

r 

7.96 ±0.10 

12.52 ±0.15 

7.11 ±0.24 

Tc [K] 

1000 

1300 

1008.68 ± 12.27 

P 

2.54 ± 0.01 

2.19 ±0.01 

2.28 ±0.01 

TL[pho] 

WcdCspe] 

22.41 

63.31 

9.86 

3.35 

3.66 

5.63 

tL[tot] 

60.13 

161.77 

65.88 

6.12 

11.56 

7.90 


Notes. Parameters without errors were fixed (not fitted) in MOD. See Sect|^for details. 


corrected ISO spectrum has its maximum flux at a 
longer wavelength than the hydrostatic model but 
the depth of the 3 jum feature is very similar to the 
model. This indicates that any difference between the 
corrected and the original ISO spectrum is caused 
by residual dust emission and not an incorrect model 
temperature. 


As there is no generally agreed value of the dust 
condensation temperature in the literature, we per¬ 
formed the fit using three different T^ values; (i) 
Tc = 1000 K, based on the theoretical equilibrium 


condensation temperature from Gail & Sedlmayr 
( 2013| l, and (ii) Tc = 1300 K (an average of the 
dust condensation temperature values from|Nowotny 


et al.|20lT and from Lobel et al.|1999| l. (iii) Finally 


we allowed all the input parameters to vary freely 
to verify the behavior of the models. The results are 
shown in Table |4] 


To check the importance of the ISO spectrum for 
the final result of the MOD fit we made some experi¬ 
ments using a coarser sampled ISO spectrum and ar¬ 
tificially increasing the error associated with the ISO 
spectrum. These tests show that the full ISO spec¬ 
trum is essential for constraining the parameters of 
the dust envelope. Nevertheless, the real uncertain¬ 
ties of the dust parameters are considerably higher 
than their formal errors. 

In Fig. we present the comparison between the 
photometry (violet circles), the spectroscopy (blue 
line) and the SED of the best fitting model. We note 
that the photometry shortward of 1 /rm cannot be 
well reproduced even if a maximum of free parame¬ 
ters is allowed. An hotter effective temperature can¬ 
not solve this discrepancy as it would be incompat¬ 


ible with the strength of the molecular features in 
the ISO spectrum. Figure shows the comparison 
between MOD visibilities (lines) and observed visi¬ 
bilities (circles). As the observed visibilities do not 
depend strongly on the wavelength we selected only 
the two wavelengths corresponding to the maximum 
and minimum visibility (Fig.[^. Comparing the spa¬ 
tial frequencies at a given visibility level, one can 
conclude that the model is roughly a factor 2 too 
small. The visibilities could be brought to agreement 
by assuming a distance of ~ 300 pc but this is un¬ 
realistically small. Therefore we have to conclude 
that the model is not extended enough or is missing 
a component. Another possibility is a difference be¬ 
tween the observed and synthetic brightness contrast 
of multiple shells. 

We will now compare the observations with the 
DMAs, which in principle would provide a more 
complete and physically consistent description of 
RU Vir. 


5. Dynamic model atmospheres 


In this section we present the analysis made by us¬ 
ing the Dynamic Model Atmospheres (from now on; 


DMA 
tra fr 
script 

) from 

Mattsson et al. (201011 and model spec- 

om |Eri 

csson et al.| (|2014ll. For a detailed de- 

ion of the modelling approach see Hofner & 

Dorfi (19971 

, Hofner (|1999|l, Hofner et al.| (|20031l. 

Matts son et 

al. (|2010|l and Eriksson et al. (|20141l). 

For applications to observations we refer to Loidl 

et al. 

19991, 

|Gautschy-Loidl et al. (2004|),|Nowotny 

et al. 

20101 

and|Nowotny et al. 1(201 1 1 . Those mod- 


els are the result of solving the system of equa- 
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Wavelength [/xm] 


Fig. 7: Comparison of observed photometry 
(Table |2.2| l-Hspectroscopy with the best fitting 
synthetic spectrum based on hydrostatic model 
atmospheres and a dusty envelope (see Table -H 
Sect.lS.ll). 



q [orcsec '] 


Fig. 8: Comparison of MOD visibilities (lines) com¬ 
pared with the VLTI/MIDI observations of RU Vir 
(circles), at 9.8 (blue) and 11.4 jim (violet). 


tions for hydrodynamics, frequency-dependent and 
spherically symmetric radiative transfer, plus a set 
of equations that describe the time-dependent dust 
grains formation, growth, and evaporation. The dy¬ 
namic model start from an initial hydrostatic struc¬ 
ture. The stellar pulsation is introduced by a ’’pis¬ 
ton”, i.e. a variable inner boundary below the stellar 
photosphere, while the dust formation (only amor¬ 
phous carbon) is evaluated by the ’’method of mo¬ 
ments” (Gauger et al. 1990[ Gail & Sedlmayr|1988| l. 


The main parameters that characterize the DMA 
are; luminosity L, effective temperature Teff (corre¬ 
sponding Rosseland radius, defined by the distance, 
from the center of the star, to the layer at which the 
Rosseland optical depth equals unity), [Fe/H], C/O, 
as well as for the dynamical aspects period P, piston 
velocity amplitude Au and the parameter f; used to 
adjust the luminosity amplitude of the model. The re¬ 
sulting proprieties of the hydrodynamic calculations 
are the mean degree of condensation, wind velocity 
and the mass-loss rate. Each model provides a set of 
’’time-steps”, representing the different phases of the 
stellar pulsation. 


We computed synthetic spectra using the COMA 


code and the associated radiative transfer (Aringer 
2000[ [Aringer et al. 2009) 1. The abundances of all 


the relevant atomic, molecular and dust species were 
computed, starting from temperature-density struc¬ 
ture, and considering the equilibrium for ionization 
and molecule formation. The continuous gas opac¬ 
ity as well as the intensities of atomic and molecu¬ 
lar spectral lines are subsequently calculated assum¬ 
ing LTE. The corresponding data are consistent with 
the ones used for the construction of the models and 
are listed injCristallo et al.|(|2007|l and [Aringer et al.j 
( |2009| ). 


SiC is added a posteriori with COMA by di 
viding the condensed material from the model into 
90% amorphous carbon using data from ([Rouleau 


|& Martin 


(Pegourie 


1991|l and 10% silicon carbide based on 
1988) . 


In order to be consistent with the model spectra 
from ( [Eriksson et al.|20f4) , we have treated all grain 
opacities in small particle limi0(SPL). The temper¬ 
ature of the SiC particles was assumed to be equal to 
the one of amorphous carbon. This is justified, since 
the overall distribution of the absorption is - except 
for the SiC feature around 11.4 jjm - quite similar 
for both species. As a consequence, the addition of 
SiC would also not cause significant changes of the 
thermal model structure. It should be noted that the 
effects of scattering are not included, since the SPL 
is adopted. 


^ An inconsistent treatment of grain opacities causes 
larger errors in the results than neglecting the grain sizes. 
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Fig. 9; Reduced vs. average outflow velocity at 
the outer boundary and vs. mass loss. Different sym¬ 
bols and colors are used for different piston veloci¬ 
ties: 2, 4, 6. 


5.1. The DMA fitting procedure 

Given the difficulties of MOD to reproduce the 
photometry, and the fact that the ISO/SWS spec¬ 
trum provides the strongest constraints for gas and 
dust, we decided to first base the comparison with 
DMA only on the ISO/SWS spectrum. Thus we 
performed a minimization in the wavelength 
range: 2.9 - 25 pm between spectroscopic observa¬ 
tion (ISO/SWS spectrum) and each of the 540 mod¬ 
els of the grid (with a total of ~ 140000 time steps). 
Based on the;^f^ we obtain a best fitting time-step for 
each model. 

In Fig.j^thejif^ vs. average outflow velocity at the 
outer boundary and vs. mass loss are shown. In Table 
l^we list the parameters of the best three time steps of 
the whole grid and the correspondingThese time 
steps cover to the 68% confidence region of;^^^. The 
relevant SEDs, without SiC added, are shown in the 
left panel of Fig.[^ The right panel shows in gray all 
the time-steps of the model that belongs to the best 
fitting time-step (x^ - 9.79). In pink this best fitting 
time-step is marked, now with the artificial inclusion 
of SiC (see also inset in Fig.[T0|). 

There is a good agreement between models and 
observations in the spectral range of the ISO spec¬ 
trum. Even based on this broad range of ~ 140000 
time-steps, it is not possible to reproduce the pho¬ 
tometry shortward of Ipm. Possible explanations for 
this disagreement are discussed in Sect.j^ 


5.2. MIDI data vs. dynamic models 


Eollowing |Paladini et al. ( ^09| l we calculated the in¬ 
tensity and visibility profiles for all the time-steps of 
the best fitting DMA derived in Sect. 5.1 We com¬ 
bined all the interferometric data (even though they 
are taken at different epochs, see Sect. |2.3| l, and we 
performed a test to obtain the best-fitting visibil¬ 
ity profile. We obtained a new best fitting time-step 
for the visibilities observations of RU Vir, different 
from the one fitting best the ISO spectrum 

The time-step best fitting the visibility is shown 
in Pig. 11 The upper panels illustrate the behavior 
of the intensity profiles at two wavelengths (8.5 and 
11.4 pm), while the lower panels show the visibil¬ 
ity vs. spatial frequencies. The phase of the best fit¬ 
ting time step is 0.24, i.e. slightly different from that 
one best fitting the SED (see Sect. |^. The best fit¬ 
ting time-step has a Rosseland radius of 2.561 mas 
(500.9 Rq) and a luminosity of 3982.8 Lq. 

The intensity are characterized by a first lobe 
(the photosphere with the hotter molecular gas) fol¬ 
lowed by two shells. The first shell is brighter at 
11.4 pm, where SiC contributes. Moreover, while at 
short wavelengths (lower left panel) the visibility can 
reproduce the data, at long wavelengths (lower right 
panel) the model is off. 

The size of the model stays the same (~ 20 AU), 
but the flux ratio between the different components 
is changing. 

The discrepancy between the model and the data 
is also observed when plotting the wavelength de¬ 
pendent visibility for the three baseline configura¬ 
tions (Pig.fT^. While at short wavelengths (8-9 pm) 
the synthetic visibilities match the level of the obser¬ 
vations, at longer wavelengths there is a slope in the 
model that is not observed in the data. The differ¬ 
ences between model and observations are more and 
more pronounced the shorter the baselines is. 


6. Discussion 

While all our attempts to reproduce the SED (ISO 
spectrum + photometry) and interferometric MIDI 
data by different models showed some agreement, 
three major discrepancies remain: (1) No model can 
describe both the optical part and the infrared part of 
the SED. This may be interpreted either as an excess 
in the optical or beyond 2 pm. (2) The model visibili¬ 
ties show higher values and/or a different behavior as 
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1 10 1 10 
Wavelength [/xm] Wavelength [/^m] 


Fig. 10: Left: Observational data (ISO spectrum in black line and photometry in violet circles), compared 
with synthetic spectra of the three best fitting time-steps (belonging to different DMA). Right: Best fitting 
time-step (pink) compared with ISO (black). The gray lines are the various other phases (time-steps) of the 
same model. 


Table 5: Three best fitting time-st eps d ynamic models from the whole grid from Eriksson et al. (20141 
resulting from comparison in Sect. 5.2 Listed are the corresponding values of the^jpToFlhe~best fitting 
time-steps, the parameters, and the phase (p. 


Teff 

[K] 

[To] 

M 

[Wol 

log g 

C/0 

AUp 

/l 

M 

[Mo/yr] 


X red[0.4-25] 

2800 

4.00 

1.5 

-0.66 

2.38 

6 

2 

3.68E-06 

6.50 

9.79 

2800 

4.00 

1 

-0.84 

2.38 

6 

2 

7.56E-06 

0.53 

11.44 

2600 

4.00 

1 

-0.96 

2.38 

6 

2 

1.44E-05 

3.40 

12.72 


a function of wavelength than the MIDI data. (3) The 
DMA predict too low fluxes beyond 14 ^m. 

Since (1) and (2) apply to all model approaches 
we believe that this is mostly caused by properties 
of the object which can neither be described by the 


rather flexible MOD approach nor by the physically 
consistent DMA (see Sect. |6.2| and |6.3| l. In the fol¬ 
lowing we will compare our results with previous 
work and discuss possible reasons for the discrep¬ 
ancies. 
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X = 8.5 [/xm] 



X = 1 1,4 [/xm] 





B[m] B[m] 


Fig. 11; Comparison of interferometric observational data for RU Vir with the modelling results based on 
the DMA best-fitting the visibilities. Up: intensity profile at two different wavelengths: 8.5 jum and 11.4 /rm. 
Down: visibility vs. baseline; the black line shows the dynamic model, the colored symbols illustrate the 
MIDI measurements at different baselines configurations. 


6. 1. Stellar parameters vs. values from the 
literature 


As described in Sect.|^and|^ stellar parameters were 
derived from the fit of synthetic spectra based on 
COMARCS hydrostatic models -i- MOD to the ob¬ 
servations (see Fig. 0, and from the study of the 
SED and MIDI data with the DMAs. 


The best fitting optical depth value resulting 
from this work is T 0.55 = 7.96 + 0.10 at d = 0.55 pm. 
Lorenz-Martins et al. (|2001|l found r = 2.5 at 1 pm. 


Ivezic & Elitzur (1995i indicate an optical depth 
value at the wavelength A - 2.2 pm of T 2.2 = 0.47. 


Assuming oc 1 /d and scaling the above literature 
results to 0.55 p we arrive at 4.5 and 1.9, respec¬ 
tively. Our value is significantly higher but our result 
is based on a broader and much better sampled SED 
which should explain the differences. 


The Rosseland radius of the best fitting time 
step(s) of ~ 2.6 mas is close to the ones derived 
from the diameter {V - K) relations of [van Belle] 
( 1999| ) and van Belle et al. ( 2013| l, namely: 2.7 mas 
and 2.81 mas, respectively. The intensity distribu¬ 
tion of the DMA (Eig. [nj has a radius of ~20 mas, 
which is comparable to the half width half maxi¬ 
mum (HWHM) resulting from the GEM-EIND fit. 
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D0-G1 




Wavelength [,u.m] 
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Fig. 12: Wavelength dependent visibilities in the MIDI range. The different panels show the three different 
baseline configurations of our observations. Models are plotted in full line, observation in dashed lines, at 
different projected baselines (see color legend). 


If we adopt these numbers for the photospheric (stel¬ 
lar) and outer envelope radius, we can conclude that 
we are observing the stratification of the star out to 
~ 6 - 7 R*. Moreover, the agreement of the GEM- 
FIND radius with the DMA radius supports the fact 
that the discrepancy between the MIDI data and the 
DMA are due to different flux ratios between the var¬ 
ious shells in the photosphere and envelope, rather 
than a global difference in size. 


The effective temperature is a parameter that 
distinguishes the different DMAs. Literature val¬ 
ues for RU Vir are in the range: [1945 - 2200] K 
([Bergeat & Chevallier 2005 [Lorenz-Martins et al. 


200 1|). These values agree with the one derived from 


the Rosseland radius of the time step fitting the SED 
(Teff = 2050 K), while the value obtained for 

the MIDI fit is slightly higher - 2545 K). 
Both temperatures are lower than the ones derived 
from fitting a static model to the C 2 H 2 feature in the 
[2.9,3.6] fim region (3000 K) and of the static ini¬ 
tial model of the best fitting DMAs (2800 K). As 
elaborated on in [Nowotny et al.| ( |2005 1 dynamic ef¬ 
fects can cause the atmospheric structure of DMAs 
to be fundamentally different from any static model, 
including the static initial model of the DMAs. 
Therefore, a T^// derived from fitting a static model 
atmosphere to observations is only very indirectly 
linked to the T^ff of the static initial model of the 
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best fitting DMA, thus a comparison between the 
two is somewhat misleading in general. In the case 
of the strongly pulsating atmosphere of RU Vir, and 
especially for the [2.9,3.6] yum region we are fitting, 
this is particularly relevant ( |Loidl et al.||1999| . We 
also note that DMAs with T^ff - 3000 K often do 
not produce the outflow required to reproduce the 
SED of RU Vir (e.g. [Eriksson et al.|2014[ ). Given the 
highly non-static atmospheric structure of the DMA, 
there is no straightforward relation between these 
different temperatures but the results indicate that the 
[2.9,3.6] yum region is a reasonable indicator for the 
(static) effective temperature. 

Einally, the mass of the best fitting dynamic 
model is 1.5 Mq. This value is in agreement with 
the 1.48 Mq obtained from the period-mass-radius 
relation of jVassiliadis & Wood| ( |1993] ). 


6.2. The shape of the SED 


Given the success of both models in reproducing 
the IR-part of the SED we consider an optical ex¬ 
cess flux more likely than an IR excess and will dis¬ 
cuss this option first. The flux difference between 
observed photometry and the models (spectroscopy) 
in the V band is at least a factor of ~ 10 in the 

SED. The quality of those used BVRI data appears to 
be good as judged from later measurements, as e.g. 
[Alfonso-Garzon et al.j ( |2012| ). This author shows a 
lightcurve which, at the phase of the ISO spectrum, 
has V ^ 11.6. This value is even brighter than the 
one used in the SED (see Table [TJ; moreover, the 
amplitude in the same band is 2 mag, and this can 
still not justify our excess in V. 

A much hotter temperature of the Mira photo¬ 
sphere can be ruled out as origin of the excess on 
the basis of the strength of the molecular features 
around 3 pm. Assuming a contribution from a com¬ 
panion would thus be the most obvious possibil¬ 
ity. However, the optical spectrum (i.e. [Barnbauin 


19921 shows no signatures from another star and no 


emission lines as e.g. typical for a symbiotic system 


( jWelty & Wade|1995|) RU Vir has a GAL EX NUV- 
flux but no FUV-flux ( jBianchi et al.|2014] l. Based on 
this flux and the predictions for white dwarf models 
( Bianchi et al.|20lT| l, one can exclude the presence 
of a white dwarf. One may also use the V and I mag¬ 
nitudes to assess the properties of a possible com¬ 
panion. However, after correcting V and I for the 
contribution from the dusty envelope as fitted with 


MOD, the resulting (V - /) and My cannot be recon¬ 
ciled with any main sequence companion. Therefore 
a stellar companion cannot explain the SED shape. 


The presence of a sub-stellar companion that is 
not seen in the optical spectrum may however imply 
the presence of a dusty disk which would produce 
excess emission in the IR. A star like V Hya would 
be a possible template ( |Sahai et al.|2008| l. The above 
mentioned GALEX flux is, however, a few orders of 
magnitude lower than what is expected for such an 
evolved star with an accretion disk. As evident from 
the MOD fit, the optical and near IR part of the SED 
do not show the amount of reddening as expected 
from the dust emission seen at longer wavelengths. 
Therefore any disk must not obscure the stellar pho¬ 
tosphere, i.e. it should be seen almost face on. This 
would also be the only disk orientation compatible 
with the lack of evidence for asymmetries from the 
MIDI data. We thus consider a disk as a very unlikely 
explanation. 


One aspect that has not been discussed so far is 
the fact that the optical photometry and the ISO spec¬ 
trum were obtained many different pulsation cycles 
apart. While we have interpolated the optical and IR 
photometry to the pulsation phase of the ISO spec¬ 
trum this cannot correct for any cycle-to-cycle differ¬ 
ences. As mentioned in Sect.|^ RU Vir has a long 
secondary period in the optical which could be an in¬ 
dication for such inter-cycle variations. Interestingly, 
also our best fitting DMA has notable inter-cycle 
variations. However, the optical and ISO data were 
obtained at similar phases of the long-term optical 
changes. Therefore, to explain the full SED with our 
best fitting DMA model, inter-cycle variations and a 
rather large phase uncertainty within a cycle have to 
be combined. 


It is known, that there are differences between 
visual phases and DMA model phases and that 
these differences as well as the effects of cycle- 
to-cycle variations increase towards shorter wave¬ 
lengths ( [Nowotny et al.|20TT I. Therefore intra-cycle 
variations together with phase assignment uncertain¬ 
ties might improve the agreement between the ob¬ 
served and the DMA SEDs but are probably not 
enough. Obtaining a broad SED within one pulsa¬ 
tion cycle would be the only way to verify this. 
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6.3. The visibilities 

While the shape of the MOD visibilities are simi¬ 
lar to the observations, the systematically lower ob¬ 
served level cannot be achieved by assuming a rea¬ 
sonable smaller distance (Sect.|^. Reducing the dust 
condensation temperature can lower the visibilities 
but produces the wrong shape. Therefore we believe 
that a simple two-component MOD model (star plus 
envelope with a single power law for the density) 
cannot reproduce the visibilities. 

With regard to the slope of the DMA visibilities 
we identify two possible explanations for the dis¬ 
crepancies observed. 

(i) The first one is model related. The models 
with high mass-loss rates produce a shell-like gas 
(and dust) component that is not observed in the data, 
i.e. the actual density distribution is smoother. This 
explanation is supported by the fact that models that 
do not produce a wind also do not exhibit the kind 
of slope seen in the best fit DMA. An example of the 
visibility of this kind of models is shown in [Sacuto 
|et al.| ( |201 l| l. A comprehensive comparison of the 
dynamic model grid with interferometric observa¬ 
tions of dusty and dust-free objects will be able to 
exclude or confirm this hypothesis. 

(ii) Besides this explanation, another possibility 
is that the environment around the star is clumpy. 
This might also explain the lack of reddening men¬ 
tioned in Sect. |6.2| Because of the presence of 
clumps, the shells appear fainter than what is pre¬ 
dicted by the model. On the other hand we do not 
observe any signal of clumps (i.e. departure from 
asymmetries) in the differential phases. Future ob¬ 
servations with the second generation VLTI instru¬ 
ment MATISSE, will provide differential phases but 
also closure phases, that will be more sensitive to 
small asymmetries. Such kind of observations will 
be able to confirm the presence of clumps. 

The differences between synthetic visibilities 
based on the DMA and MIDI visibilities around 
the SiC feature are most pronounced for the longest 
baselines, i.e. closest to the star. Since the formation 
process of SiC is not really known, the calculations 
assume that the SiC abundance scales with the amor¬ 
phous carbon abundance as provided by the DMA. 
The above differences indicate that SiC could form 
further in and/or in a smaller amount but this has to 
be confirmed by comparisons based on a larger sam¬ 
ple of stars. 


The fact that the model phases best fitting the 
SED and the visibilities are different is probably 
caused by the above mentioned cycle-to-cycle vari¬ 
ations of the best DMA model. In the MIDI-range 
these variations are comparable to the intra-cycle 
variations, i.e. a certain pulsation phase within one 
cycle may give similar visibility profiles as a differ¬ 
ent phase in another cycle. 

6.4. The flux beyond 14 fim 

The long wavelength region of the SED can be well 
reproduced with MOD, while the DMA always pre¬ 
dict a too steep decline in flux towards the longest 
wavelengths. Tests showed that an extension of the 
DMA further out cannot account for this difference. 
Assuming a smoother density distribution of the 
DMA, as already argued from the visibilities and as 
present in the MOD model, should increase the dust 
emission at longer wavelengths. A more extensive 
comparison between DMA and a larger sample of 
C-rich Miras (Rau et al., in prep.) is needed to check 
whether the difference beyond 14 /urn is found only 
in RU Vir or is a general characteristic of DMAs. 

A higher mass loss in the past might be an ex¬ 
planation although the shell parameters from MOD 
are close to a stationary wind solution. If the mass 
loss has even stopped in the recent past, this might 
account also for the overall differences in the shape 
of the SED as a very small mass loss now should 
shift the emission from circumstellar dust to longer 
wavelengths. However, besides the good fit with a 
stationary wind, a higher mass loss in the past would 
affect the visual light curve. But the average visual 
magnitude and the period have been stable over the 
last ~ 100 years and therefore such a scenario is also 
not very likely. 

7. Conclusions 

As demonstrated in this study, the joint use of pho¬ 
tometry, spectroscopy and interferometry and a com¬ 
parison with models is essential to achieve a full un¬ 
derstanding of the atmospheres of AGB stars. 

In this work we presented a study of the atmo¬ 
sphere of the Mira C-rich star RU Vir. We com¬ 
bined spectroscopic photometric and interferomet¬ 
ric observations, and hydrostatic and dynamic model 
atmospheres as well as simple geometrical models. 
Until now, the only studies with a similar approach 
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as in the present study are the one by Sacuto et al. 
( |20 TT] i on R Scl and RT Vir [Sacuto et al.| P013| l 
on RT Vir. While R Scl is a Semi-regular pulsat¬ 
ing star, and RT Vir is an M-type star, RU Vir is the 
hrst carbon-rich Mira variable for which photome¬ 
try, spectroscopy and interferometry are compared to 
DMAs. Studying C-rich stars is particularly impor¬ 
tant since DMAs are more advanced for the C-rich 
stars than for the O-rich stars, because the dust for¬ 
mation process is considered thought to be simpler 
and better understood in the C-rich case (jNowotnyj 

|etal.|2010l|20TTll. 

The HWHM derived by the GEM-FIND htting 
shows overall agreement with the DMA size in the 
V-band (see Sect. EB- Furthermore the Rosseland 
radius corresponding to the time-step best htting the 
interferometric data is in agreement with the diam¬ 
eter derived via the (V - K) relation by van Belle 
( |1999| l; [van Belle et al.j ( |2013] l. The htted effective 
temperature of the best-htting DMA time steps is 
in the range of previous observational estimates but 
much lower than the temperature of the hydrostatic 
model of Sect. 14. H and of the initial model of the best 
htting DMA. 

The shape of the SED in the ISO range can be 
well reproduced both with MOD or DMA, some 
discrepancies remain shortward of 2 /am and long- 
wards of 14 fj.m. A similar situation exists with re¬ 
gard to the interferometric data, both in the shape 
and in the visibility level. Some of the differences 
might be explained by a combination of intra- and 
inter-cycle differences as the observations are spread 
over many pulsation cycles and both the star and the 
best htting DMA show inter-cycle variations. Other 
possible reasons could be a decrease in mass loss 
over the last few hundred years or a sub-stellar com¬ 
panion associated with a dusty disk. However, these 
scenarios are not considered to be very likely and 
to check them both further observations and mod¬ 
elling are necessary. Based on the current work, we 
suspect that RU Vir is a somewhat peculiar object. 
Extending the comparison between models and ob¬ 
servations to a larger sample of C-rich Mira variables 
(Rau et al., in prep.) will be necessary to clarify this, 
and will provide the general characteristics of the at¬ 
mospheres of these stars and further constraints for 
the models. 
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Appendix A: Observing iog 
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Table A.l: Journal of the MIDI observations of RU Vir. 


Target 

UT date & time 

Config. 

Sp 

[m] 

PA 

[°] 

Seeing 

[”] 

Airmass 

Mode 

<!> 

RU Vir 

2010-05-23 T01:23:12 

HO-EO 

47.2 

73 

1.14 

1.52 

SCI-PHOT 

0.11 

HD 120323 

2010-05-23 T00:32:16 


45.6 

53 

1.14 

1.12 



HD 81797 

2010-05-23 T23:05:53 


47.4 

73 

0.99 

1.06 




2010-05-23 T00:00:32 


44.6 

74 

0.70 

1.18 



HD 133216 

2010-05-23 T02:18:36 


65.0 

98 

1.14 

1.05 



RU Vir 

2010-06-04 T01:30:24 

HO-EO 

47.9 

72 

1.45 

1.17 


0.14 

HD 120323 

2010-06-04 T02:06:52 


43.8 

45 

1.43 

1.21 




2010-06-04 T00:50:11 


47.6 

63 

1.41 

1.03 




2010-06-04 T03:03:39 


45.2 

82 

1.20 

1.07 




2010-06-04 TOl: 12:54 


47.9 

67 

1.63 

1.02 




2010-06-04 T01:50:48 


47.8 

72 

1.53 

1.02 



HD 133216 

2010-06-04 T03:46:05 


47.0 

68 

1.43 

1.17 



RU Vir 

2014-04-11 T00:48:17 

DO-Gl 

67.6 

134 

0.96 

1.52 

HIGH-SENS 

0.37 

HD 120323 

2014-04-11 T01:03:12 


45.3 

112 

0.79 

1.70 




2014-04-11 T01:28:ll 


49.7 

113 

0.96 

1.52 



RU Vir 

2014-04-11 T01:15:18 


69.4 

132 

1.00 

1.59 


0.37 

HD 120323 

2014-04-11 T01:03:12 


45.3 

112 

0.79 

1.70 




2014-04-11 T01:28:ll 


49.7 

113 

0.96 

1.52 



RU Vir 

2014-04-11 T02:11:25 

HO-H 

40.5 

141 

1.02 

1.32 


0.37 

HD 120323 

2014-04-11 T02:25:17 


35.2 

134 

0.95 

1.26 



HD 81797 

2014-04-11 T02:55:16 


36.5 

163 

1.03 

1.21 




2014-04-11 T03:20:45 


36.0 

168 

1.13 

1.30 



RU Vir 

2014-04-11 T02:38:08 

HO-H 

40.0 

142 

0.98 

1.24 


0.37 

HD 120323 

2014-04-11 T01:58:26 


33.7 

133 

0.80 

1.36 




2014-04-11 T02:25:17 


35.2 

134 

0.95 

1.26 



HD 81797 

2014-04-11 T02:55:16 


36.5 

163 

1.03 

1.21 




2014-04-11 T03:20:45 


36.0 

168 

1.13 

1.30 



RU Vir 

2014-04-11 T03:07:40 

HO-H 

39.3 

143 

1.05 

1.18 


0.37 

HD 120323 

2014-04-11 T02:25:17 


35.2 

134 

0.95 

1.26 



HD 81797 

2014-04-11 T02:55:16 


36.5 

163 

1.03 

1.21 




2014-04-11 T03:20:45 


36.0 

168 

1.13 

1.30 



RU Vir 

2014-04-11 T03:58:15 

DO-Gl 

68.4 

130 

1.18 

1.14 


0.37 

HD 120323 

2014-04-11 T03:45:37 


66.7 

124 

1.32 

1.07 




2014-04-11 T04:ll:45 


68.5 

127 

1.20 

1.04 




2014-04-11 T04:39:44 


69.9 

130 

1.38 

1.02 




2014-04-11 T05:47:28 


71.5 

138 

1.12 

1.03 



RU Vir 

2014-04-11 T04:25:14 

DO-Gl 

66.2 

131 

1.28 

1.14 


0.37 

HD 120323 

2014-04-11 T03:45:37 


66.7 

124 

1.32 

1.07 




2014-04-11 T04:ll:45 


68.5 

127 

1.20 

1.04 




2014-04-11 T04:39:44 


69.9 

130 

1.38 

1.02 




2014-04-11 T05:47:28 


71.5 

138 

1.12 

1.03 



RU Vir 

2014-04-11 T04:53:53 

DO-Gl 

63.4 

133 

1.26 

1.16 


0.37 

HD 120323 

2014-04-11 T03:45:37 


66.7 

124 

1.32 

1.07 




2014-04-11 T04:11:45 


68.5 

127 

1.20 

1.04 




2014-04-11 T04:39:44 


69.9 

130 

1.38 

1.02 




2014-04-11 T05:47:28 


71.5 

138 

1.12 

1.03 



RU Vir 

2014-04-11 T05:19:21 

DO-Gl 

60.6 

135 

1.10 

1.20 


0.37 

HD 120323 

2014-04-11 T04:11:45 


68.5 

127 

1.20 

1.04 




2014-04-11 T05:47:28 


71.5 

138 

1.12 

1.03 



RU Vir 

2014-04-11 T05:59:56 

DO-Gl 

55.6 

139 

1.09 

1.30 


0.37 

HD 120323 

2014-04-11 T04:11:45 


68.5 

127 

1.20 

1.04 




2014-04-11 T05:47:28 


71.5 

138 

1.12 

1.03 
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Fig. A. 1: Comparison of the observational MIDI data with the best-fitting geometrical model (See Table [^. 
Shown are visibilities vs. wavelength at different position angles and baseline lengths (See Fig. I^for the 
complete set of visibilities). Black lines represent the MIDI observations of RU Vir. Light blue lines represent 
the fit of the best-fitting UD-tGaussian model. 














